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ABSTRACT 

Using a simple two-dimensional, zero-/? model, we explore the manner by which reconnection at a current 
sheet releases and dissipates free magnetic energy. We find that only a small fraction (3%-l 1% depending on 
current sheet size) of the energy is stored close enough to the current sheet to be dissipated abruptly by the 
reconnection process. The remaining energy, stored in the larger-scale field, is converted to kinetic energy in 
a fast magnetosonic disturbance propagating away from the reconnection site, carrying the initial current and 
generating reconnection-associated flows (inflow and outflow). Some of this reflects from the lower boundary 
(the photosphere) and refracts back to the X-point reconnection site. Most of this inward wave energy is 
reflected back again, and continues to bounce between X-point and photosphere until it is gradually dissipated, 
over many transits. This phase of the energy dissipation process is thus global and lasts far longer than the initial 
purely local phase. In the process a significant fraction of the energy (25%-60%) remains as undissipated fast 
magnetosonic waves propagating away from the reconnection site, primarily upward. This flare-generated 
wave is initiated by unbalanced Lorentz forces in the reconnection-disrupted current sheet, rather than by 
dissipation-generated pressure, as some previous models have assumed. Depending on the orientation of the 
initial current sheet the wave front is either a rarefaction, with backward directed flow, or a compression, with 
forward directed flow. 

Subject headings: MHD — shock waves — Sun: flares 



1. INTRODUCTION 

Magnetic reconnection has been frequently proposed as 
the mechanism whereby magnetic energy, stored in the so- 
lar corona, is rapidly released in a solar flare. In most cur- 
rent models, fast magnetic reconnection occurs at a current 
sheet where magnetic field lines of differing connectivity are 
brought into close enough proximity for a small-scale process, 
such as Ohmic diffusion or various kinetic effects, to forge 
new connections between them. The current sheet, by carry- 
ing a net current, also stores the magnetic energy which the re- 
connection liberates. This energy is not, however, co-located 
with the current sheet itself so the mechanism responsible for 
the reconnection electric field will not be the same one respon- 
sible for releasing or dissipating the magnetic energy. Many 
investigations have been focussed on the former, local mech- 
anism (flux transfer) while far fewer have addressed the latter, 
global mechanism (energy release). 

A concrete illustration of the above issue is provided by the 
simple, two-dimensional quadupolar coronal field shown in 
Fig. [TJ A pair of sources, P2 and N2, have emerged under- 
neath an older bipole, Pl-iVl. If not all of the new flux is able 
to reconnect with the overlying flux, the state of minimum 
magnetic energy, B(x), will contain a curren t sheet separat- 
ing new from old flux as sho wn in Fig.QJ (see Hevva rts et alj 
1977; Priest & Forbes 2000). The state of lowest possible en- 
ergy is a potential field, Bo(x), shown in Fig.QJ, in which 
additional flux Aip (grey regions) interconnects the new and 
old polarities. A current sheet, carrying current I cs , exists to 
maintain a connectivity different from the potential field, as in 
Fig. [TH . For this reason the net cu rrent I cs ~ Aip (Longcope 
l200lFlLongcope & Magarall2004 . The current sheet will be 



at the current sheet may have energetic consequences. 



reduced, or eliminated entirely, as flux is transferred across 
the sheet into domain connecting P1-N2 and P2-N1 (the 
shaded regions). It is in this way the local topological changes 




FIG. 1 . — The field from four photospheric, magnetic sources, (a) Equilib- 
rium B(x) state after emergence and incomplete reconnection of the bipole 
P2-N2. The current sheet is shown as a bold arc connected to thinner lines: 
the separatrices. (b) The potential field, Bo(x), accessed by reconnection 
transfering flux Ai/> into connection between P2-NI andPl-7V2; the new 
field lines fill the shaded regions, (c) The non-potential field, B — Bo, due 
to the current sheet. 

The free magnetic energy of the topologically constrained 
magnetic field is computed by subtracting the energy of the 
unconstrained (i.e. potential) field, 



1 



8tt 
1 

8^ 



AE M = — / \B\ z d''x - — 



8tt 



\B \ 2 d 3 x 



|B-B | 2 d 3 a; 



(1) 



where the integrals are over the entire coronal volume (z > 



2 



0)Q The final expression shows that the free energy is equiv- 
alent to the energy of the non-potential component, B(x) — 
Bo(x), due to the current sheet with homogeneous conditions 
at the lower boundary (B z — Bq z = 0). This energy will natu- 
rally decrease as the current I cs is decreased by reconnection. 
The non-potential field, shown in Fig.[U;, extends far beyond 
the sheet itself. At large distances it decreases inversely with 
distance similar to the field from a simple wire. 

It is evident from the form of the non-potential field (Fig. 
QJ;) that changes in the current sheet must propagate great dis- 
tances in order to release the stored magnetic energy. Much 
of the volume over which the free magnetic energy is stored 
is not in magnetic contact with the current sheet, so neither 
Alfven waves nor slow magnetosonic waves (nor shocks) can 
change the field directly. Yet if the current in the sheet is to 
change, that change must be reflected in the distant field, in 
accordance with Ampere's law. How this occurs is unlikely 
to depend too critically on the mechanism responsible for the 
flux transfer which must, according to all current theories, oc- 
cur on small scales. Thus it should be possible to study the 
mechanism of energy release using a arbitrary form of rapid 
flux transfer at the current sheet. 

Such an analysis was performed by Longcope & Priestl 
(2007) using a current sheet situated at a two-dimensional X- 
point, in a zero-/? plasma, whose reconnection was effected 
by a sudden increase in Ohmic resistivity. They found that the 
disruption of the current sheet launched a cylindrical current 
shell at the leading edge of a fast magnetosonic wave (FMW). 
The shell contained almost all of the current formerly carried 
by the current sheet. It therefore left in its wake (i.e. inside 
the shell) a nearly potential magnetic field. The non-potential 
magnetic energy of the initial field was converted to kinetic 
energy of the wave's flow. Remarkably, the Ohmic dissipa- 
tion responsible for initiating the wave, and thus releasing the 
stored magnetic energy, directly dissipated a rather small frac- 
tion of the stored energy. This puzzling result, that large resis- 
tivity does not result in large Ohmic losses, can be anticipated 
from the arguments above: even after it diffuses, the current 
sheet occupies a very small volume and therefore has access 
to very little magnetic energy. 

The flow field in the FMW has a quadrupolar structure fa- 
miliar in steady reconnection models: inflows along one axis 
(the vertical axis in Fig.Q]) and outflows along the other (hor- 
izontal). Previous studies of transient magnetic reconnection 
models have found fast magnetosonic ra r efaction waves to 
be th e drivers of inflow dLin & Ledll994t iHeyn & Semenovl 
1996; Nitt a et ail 120011) . These analyses, set on infinitely 
long current sheets, focussed on the inner reconnection region 
rather than the front of the wave. Since Longcope & Priestl 
d2007l) considered a sheet carrying finite current they were 
able to analyze its full global propagation and energy release. 
In particular, their model revealed that the kinetic energy of 
the inflow must be supplied by decreasing the free magnetic 
energy stored in the initial equilibrium. In their unbounded 
domain this kinetic energy is not subsequently converted to 
heat. 

The finite current sheet studied by Longcope & Priestl 
(2007) was situated in an unbounded domain, which therefore 
contained an infinite amount of free magnetic energy. The 
FMW would propagate indefinitely, converting this stored en- 
ergy to kinetic energy at a uniform rate. The wave energy 

1 A cross term involving the integral of (B — Bo) ■ Bo can be seen to 
vanish after expressing Bq = — Vx and integrating by parts. 



could therefore become arbitrarily large in comparison to the 
energy directly dissipated at the current sheet. 

In order to mak e conta ct with previous w ork b y 
[ Craig & McClvmontl d 199 11) and iHassaml dl992l) , 
Longcope & Priest (2007) briefly considered the effect 
of a boundary: a concentric cylindrical conductor. This 
reflected (perfectly) the outward-propagating FMW back 
inward. T he inward wave collapsed on t he X-po i nt, as 
described bv lCraig & McClvmontl d 19911) and lHassaml d 19921) 
permitting still more Ohmic dissipation to occur there. The 
majority of the wave's energy was, however, reflected once 
more by the X-point. The wave was therefore trapped 
between a perfectly reflecting outer conductor and an imper- 
fectly reflecting X-point. The only losses were at the X-point, 
so eventually all energy was in fact dissipated at the X-point 
through Ohmic dissipation. The characteristic time for 
dissipation depended on the round-trip transit time between 
reflectors. Owing to the exponential increase of the Afven 
speed with distance, this transit time scales logarithmically 
with the dissipation scale, and thus logarithmically with the 
resistivity. The ability of an X-point to Ohmically dissipate 
magnetic energy in logarithmic time was the most sign ificant 
result of the studies Crai g & McClvmontl ( 1199 11) and lHassaml 
dl992l) . and appears to suggest that Ohmic dissipation alone 
is capable of all magnetic energy dissipation in a flare. 

The upshot of this reasoning is that the initial disruption of 
the current sheet by Ohmic diffusion will dissipate a small 
fraction of its magnetic energy, but that repeated reflections 
can result in the complete dissipation of all stored energy at 
the X-point. To achieve the latter end, however, previous au- 
thors have assumed a conducting boundary completely sur- 
rounding the X-point. In a more realistic geometry, such as 
that of Fig. Q] there is a conducting boundary at the photo- 
sphere (z — 0) but it does not completely surround the X- 
p oint. 

iMcLaughlin & Hoodl d2006a) studied the linearized, ft — 
dynamics of a quadrupolar magnetic field with a planar lower 
boundary like Fig. [T] Rather than initiate a wave by recon- 
nection, they launched a FMW by fiat from the lower bound- 
ary. They found that a fraction of the wave (they report 40%) 
was refracted into the null, while the remainder continued to 
propagate away. This suggests that the photospheric lower 
boundary may indeed be less effective at mediating energy 
dissipation than are the concentric cylindrical conductors, al- 
though it is not clear that the same fraction would apply to 
a wave launched outward from the null rather than from the 
lo wer boundary. 

McLaugh lin & Hoodl (|2006a) omitted diffusion from their 
model but assumed its effect would b e to dissipate all the 
wave energ y refract e d into the null. [Craig & McClvmont 
d 199 lb and IHassaml d 19921) found to the contrary that the 
main effect of resistivity is to reflect the wave back out- 
ward from the null. It was only through repeated reflec- 
tion that the wave energy could be ultimately dissipated. If 
on ly 40% were reflected back at each step (as suggested 
by McLaughlin & Hood (2006a)), then the final dissipation 
would be far less than 40% initially directed toward the null. 

The lack of 100% reflection back to the null point appears 
to pose a difficu lty for the system r e aching a potential state. 
It was found by Longco pe & Priestl (120071) that flux transfer 
continues at the X-point long after the majority of the current 
sheet had disappeared. As a result Aip continues to decrease, 
below zero, bringing the system away from the potential field 
state (Alp — 0). When the wave reflected from the cylindrical 
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boundary, however, it reversed this flux transfer, bringing Aijj 
back upward. While it continued to overshoot the potential 
value, numerous reflections led to its gradual convergence to 
zero. If the initial wave is only partially reflected back to the 
null point, as the results of McLaug hlin & Hoodl d2006al) sug- 
gest, then it is unclear how Aip would ever converge to zero, 
and the system achieve a potential field. 

To more fully understand energy release and dissipation we 
seek to determine what fraction of the FMW launched from 
the null point will be reflected back to the X-point in this 
more realistic configuration. We must also know how much 
of the partially reflected wave will ultimately dissipate at the 
X-point and what fraction will reflect once more. 

This analysis, performed below, shows that the lowest fre- 
quency components of the wave are reflected almost com- 
pletely by the photospheric boundary. This results in the cur- 
rent at the X-point being entirely eliminated, and At/j — > 0, 
after numerous reflections. It also results in the direct dis- 
sipation of a significant fraction of the free energy (between 
40% and 75% in one example we consider). This does not oc- 
cur immediately through the local reconnection mechanism, 
Ohmic diffusion in our case, but rather it requires repeated 
reflections from the photosphere and thus takes many transit 
times to achieve. The remainder of the free energy (25% - 
60%) is emitted by the reconnection as FMW at higher fre- 
quencies. These waves form a number of pulses propagating 
primarily vertically upward, away from the photosphere. The 
flow in these wave forms is either a rarefaction or a compres- 
sion depending on the orientation of the initial current sheet 
(horizontal or vertical, respectively). Such fast magnetosonic 
disturbances have long been known to accompany flares, but 
until now no quantitative reconnection model has predicted 
what fraction of the free magnetic energy they accounted for. 

We present the model calculation for a two-dimensional, 
quadrupolar field, like the one in Fig.[TJ In the next section we 
specify the geometry of the model field, and quantify the free 
magnetic energy stored in advance of reconnection. In the 
following section we analyze the dynamical behavior in the 
vicinity of the X-point, and describe how this is matched to 
the external field, including the photospheric boundary. Sec- 
tion [4] presents numerical solutions to the external response, 
including the emission of a FMW upward. The numerical so- 
lutions cannot be continued for the many FMW-transit times 
of the full solution. Instead we use them to characterize the 
reflection from the photosphere and then use the reflection 
coefficient to synthesize a solution for long times in £|5] This 
long time solution is used to quantify the ultimate fate of the 
free energy. 

2. MODEL OF RECONNECTION IN A QUADRUPOLAR FIELD 

2.1. The quadrupolar equilibrium 

We begin with a potential magnetic field created by four 
photospheric sources, as depicted in Fig. [TJi. The field is ex- 
pressed in terms of a flux function B(y, z) = VA x x, whose 
potential version is 
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where source j is located at (yj, —d) and has flux (per ignor- 
able length) ipj . Singular line sources are located a distance 
d below the photospheric boundary (z = 0) so that the field 
is non-singular within the model (z > 0). The sum runs over 



positive sources j = 1, 2 and negative sources j = —1, —2. 
We create reflectional symmetry about the z axis by taking 
ip-l = —ipi and -0_2 = — ip2 as we ll as V-x = —yi and 
j/_2 = —yi- Except for scaling, the quadrupole is described 
by the two dimensionless parameters ipi/ip2 and 2/2/2/1; we 
take the latter to fall in the range (—1, 0) so that j — —2 and 
+2 create the inner bipole whose orientation is opposite to the 
outer bipole. 

The potential field generated by (fJI will vanish at a single 
coronal X-point located at (0, z x ), where 
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provided this is a real positive number. There are two condi- 
tions for the expression to be real. First, that the denominator 
inside the radical must remain positive, I2/2I/2/1 < V'i/^. 
which is equivalent to requiring the overlying field to be di- 
rected opposite to the lower bipole. The second is that the 
numerator be positive, I2/2I/2/1 < ips/ipu m order that the the 
field direction not reverse along the line between P2 and N2. 
(When this is violated two null points are located symmetri- 
cally along the line between P2 and N2 — at z = —d.) Once 
these conditions are met it is still necessary that z x > in 
order that the null point is actually above the lower boundary 
at z = 0. 

In the immediate vicinity of the null point the potential field 
flux function, eq. @, is approximated by its Taylor expansion 



A (y,z) ~ A (0,z x ) + ±B' [y 2 -{z-z x ) 2 
where the magnetic shear at the X-point is 
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This i s the same structure analyzed by Longcope & Priestl 
(l2007h . and we will use those results over an inner region, 
extending to a small distance from the null point. Within 
this region we assume the expansion © to be a good approx- 
imation to the full field, eq. y). 

We add to the potential field, Bo, the non-potential field 
from a current sheet carrying current 
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The full field, B(x), is a force-free equilibrium with a singular 
current, similar to that in Fig. QJ. The current will be taken 
small enough that the sheet is a straight line of half-length 

A = 2y/\I cs \/B' , centered at (0, z x ). If 7 CB > the sheet 
will be horizontal, similar to Fig. [T^; if 7 CS < it will be 
vertical. We will assume the current is small enough that A <C 

r<. 

In order to preserve the distribution of photospheric flux 
an image current sheet, with opposite sign, will be added at 
z = —z x . The com bination will yie ld a flux function value 
on the current sheet (Longcope 2001) 



A cs ~A o {0,z x ) + J cs ln(16ez2/A 2 ) 
= A o (0,z x ) + 7 cs ln(4z 2 x eB /|/ cs |) 



(7) 



where e is the base of the natural logarithm. The result of 
the current sheet is thus to introduce an additional flux At/j = 
I cs ln(4z 2 ,eB /|/ cs |) connecting P2 to N2. The smallness of 
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the current sheet means that outside the inner region radius 
(fi) the full field may be approximated 



A(y,z)~A (y,z) 
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y 2 + (z + Zy) 2 
A (y,z) + 2/ cs ln(2z x /r) , 
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where r is the distance from (0, z x ). 

The free energy of the magnetic field can be computed by 
integrating the electroma gnetic work requ ired to ramp the cur- 
rent up to its final value (Longcope 200l]) 
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The amount of this free energy outside the inner region is 



AE (out) = ^_ 

8tt 



- / \VA- VA \ 2 d 2 x 

>7T J r>n 



~±/ 2 s l n (4z 2 /r 2 ) . (10) 

It is due to the lower boundary, and the image current, that 
these e xpressions are finite, while those of Longcope & Priestl 
(2007) were infinite. 

2.2. Linear dynamics 

Having assumed a small current sheet we can consider its 
contribution to be a s mall perturbation, A^ , to the potential 
quadrupolar field (see Longcope & Priest 2007, for extensive 
discussion of this approximation). This evolves, along with 
a linear velocity field, vi, according to the momentum and 
resistive induction equations of resistive MHD 



^i = _Z^£v 2 A 1 
dt 4:np Q 
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= -Vi- VA + r^ 2 A 1 , 
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where r\ is the Ohmic resistivity and po the plasma's initial 
mass density. The pressure force has been dropped under the 
assumption t hat plasma is very low. 

Following Longco pe & Priestl (120071) . we introduce an in- 
ductive electric field variable, U± = vi • VAo = — x ■ (vi x 
Bo), to obtain the pair of scalar equations 
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where va = \VAo\/y/4wpo is the Alfven speed in the poten- 
tial field. For simplicity we henceforth take the initial mass 
density, po, to be uniform. The system is solved beginning 
with Ai(y, z, 0) generating the current sheet. The resistivity 
is then "turned on" to some fixed value at t = 0. The current 
sh eet then diffuses ini t iating the dynamical evolution studied 
by lLongcope & Priestl (120071) . 

Within the inner region the potential magnetic field in- 
creases linearly with distance from its X-point, 
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r = vV + {z- z x ) 2 < r z . (15) 



This means the Alfven speed increases linearly as well, va = 
uja t, where lua — B' Q / 'y/inpo is a frequency characteristic 
of the null point. This frequency defines a resistive length 
scale £ v — y/rj/tJA, inside which diffusion is a dominant ef- 
fect. We propose that the enhanced resistivity turned on at 
t = be large enough that £ v ^> A. This will assure that 
the current sheet is entirel y disrupted by the diffusive process 
(Long cope & Pri est 2007). This is the way we model simply 
a mode of fast magnetic reconnection in order to focus on the 
external response. 

While the diffusive length is much larger than the current 
sheet, we take it to be small compared to the potential field 
X-point. This permits us to define the inner region such that 
r% 3> £ v , and the resistive term may be dropped from eq. ( TBI 
outside of it ( r > rj). We will then be a ble to match the in- 
ner solution of Long cope & Priesfl d2007l) to the non-diffusive 
solution of the outer region. We will first consider, in the next 
section, solutions in the interface region r ~ ri at which the 
matching must occur. We then produce the complete solution 
of the inner region in a form suitable for matching. 

3. DYNAMICS OF THE INNER REGION 

3.1. The interface region 

In the vicinity of the interface, r = ri, we may neglect 
the resistive term (r, 3> £ v ) and use the X-point expansion 
of the Alfven speed. We take advantage of the approximate 
cylindrical symmetry to expand 
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where the sum extends over both positive and negative values 
of m. The modal coefficients then satisfy the equations 
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The axisymmetric equations (m = 0) are satisfied by arbitrary 
wave forms 
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propagating outward (W+) or inward (W-). Because the 
wave speed increases linearly with radius, a particular phase 
point in the wave accelerates exponentially as it travels: r ~ 

Higher mode numbers (|m| > 1) satisfy a Klein-Gordon 
equation in In r, with a cut-off frequency, |to|o;a. Frequencies 
below the cut-off are evanescent at the interface. We demon- 
strate below that higher-m perturbations of low frequency in- 
cident on the X-point will be reflected from the interface. 

The initial magnetic field in the interface region is given by 
the second term of eq. (|S), for which rdA\ /dr = —2I CS . This 
is an axisymmetric field in equilibrium ( Ui = 0) consisting of 
counter-propagating waves of the form 



W-{x) 



W+(X) = 2U A I C : 



(21) 



It seems that in the absence of resistivity the X-point and cur- 
rent sheet function as a perfect reflector, at least at zero fre- 
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quency. It reflects an incoming wave with a sign flip in order 
to maintain the current-sheet equilibrium. We see below that 
the onset of resistive diffusion changes the nature of this re- 
flection, thereby eliminating the current sheet. 

Integrating the Poynting flux around the entire the interface 
at fj yields the power carried outward by the waves 



The axisymmetric contribution to this power 
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is the difference in energy carried in and out across the inter- 
face. 

3.2. The inner region: current disruption 

Longcope & Priestl (12007b performed a thorough analysis 
of the dynamics following the diffusive disruption of the cur- 
rent sheet; pertinent results are summarized in Fig. [2] The 
external effect of the current disruption is an axisymmetric 
FMW propagating outward. The axisymmetric electric field, 
tj( \ at a fixed position (Fig. [2^) rises from zero as the FMW 
reaches it. They found the form of the wave to beQ 
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where to is the time required for the wave to be initiated and 
reach the nominal position r = 1. This semi-empirical fit, 
plotted as a dashed curve over the actual wave form in Fig. 
[2J3, is a reasonably good match. 
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FIG. 2. — External response to current sheet dissipation after 
Longcope & Priest 12007). (a) The axisymmetric electric field, U^°\ scaled 
to oja^cs, at r = 12£ v (solid) and r = I5OA7 (dashed) for comparison, 
(b) The outgoing wav e, W + at r = 12£ v (solid), scaled to uip^I CB , and the 
fit (dashed) from eq. i24\ using the value of to shown by the arrow. (This 
is the only free parameter in the fit.) (c) The axisymmetric contribution to 
the Poynting flux, scaled to u^I^ B , at r = I2£ v (solid). Its time-integral 
(dashed), the total energy transferred, is scaled to and plotted against the 
axis on the right. 

2 There is a typographical error in the text following eq. (30) in 
Longcope & Priest 12007). It should have read rrj = V2£ v . Equation i24\ 
corrects that error. 



Their analysis assumed an infinite exterior domain from 
which a steady incoming wave, W- = 2u>aIcs, persisted 
indefinitely. The outgoing wave in eq. d24t . switches from 
being negative, opposite to the incoming wave, to positive, 
W+ — > 2ljaI cs , matching the incoming wave (see Fig. Wp)- 
This represents a change from a static magnetic field with no 
velocity or electric field (an equilibrium with U\ = 0) to a 
steady electric field, XJ\ = 2uj^I cs (see Fig.|2^), and vanish- 
ing magnetic field perturbation (i.e. the potential field). Evi- 
dently the change in resistivity (i.e. its turn-on) produces this 
change in the X-point reflection properties. 

The Poynting flux becomes negative as the FMW passes a 
given radius (see Fig. [2};). It can be seen from eq. ( 1231 that it 
reaches a peak value, = —luaI 2 s /2, at the instant W+ = 0. 
After that W+ —> W- causing the Poynting flux to decrease 
asymptotically as t -1 . The same inward Poynting pulse is ob- 
served, with some travel delay, at all distances from the cur- 
rent sheet independent of the value of resistivity responsible 
for reconnection. 

There is a final aspect of this wave solution worthy of 
note. As the solution approaches the steady electric field, 
Ui — > 2uj\I cs , eq. ( fl9b implies the flux function will become 
uniform, but with a steadily decreasing value 
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This steady flux transfer, dictated by the external solution, 
is accommodated by a persistent current density confined to 
an ever-shrinking neighborhood of the X-point well inside 
the diffusive radius £ v . The central current density remains 
fixed, so the total current decreases as t -1 , the square of the 
shrinking diameter. This current, amidst a steady electric field 
U ~ 2wA-fcs, produces the asymptotic Poynting flux seen in 
Fig. |2};. The FMW imposes an inflow and the residual internal 
solution represents its accommodation at small scales. 

3.3. The inner region: response to outer dynamics 

The complete dynamics of the inner region, including resis- 
tive diffusion, is best treated after Fourier transforming 
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The Fourier coefficient must then satisfy the equation 
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iHassaml ([1992) found an analytic solution for the case m — 
which was regular as r —> 0. This solution can be generalized 
to arbitrary m, 



r m F[ a,a*:m + l; 



MAT 

where F is the hypergeometric 
( Abramowitz & Stegun 1964) and 



1 r 
2 



. + i\J (uj/uja) 2 — m 2 



lVJ r 



(28) 
function 

(29) 



For frequencies above cut-off, uS 1 > m 2 u! A , the factor -crj m 
is real and the hypergeometric function can be expanded well 
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outside the diffusion region, r ^> £ rp 



' W A 

XI 



-m/2 



.WAT 



,/2 



»/2 



(30) 



where X m = m\T(iwm) I 'aT 2 (a). The factors 



= exp 



±im m Infr/ij]) ± \iw m ln(zwA/w) 



(31) 



can be identified with outward (upper sign) and inward (lower 
sign) propagating components, assuming that w m is real and 
positive. An analytical continuation through complex fre- 
quencies from one branch of the expression onto the other 
would require an exchange between what have been here des- 
ignated as "inward" and "outward" components, thereby in- 
verting the expression for the reflection coefficient defined be- 
low. 

The ratio of the terms corresponding to outward and inward 
waves in eq. ( f3Qb gives a complex reflection coefficient 

X,, 



R {m \uj) = -^- exp[im m \n(iui A /oj)} 



X„ 
A 7 * 



(32) 



where a minus sign is introduced to represent the reflection of 
Ui rather than A\. The magnitude of the coefficient represents 
the ratio of amplitudes between the incoming and the reflected 

wave, 



\R (m) (u)\- 



;exp 
= exp 



;VJ ri 



(33) 



Resistivity at the X-point damps all waves above the cut-off 
frequency: \R\ < 1. At frequencies approaching the cut-off 
the X-point becomes a more perfect reflector; at frequencies 
below the cut-off it is a perfect reflector. 

The current sheet provides an initial condition 
A-\ (r, (j), 0) containing all even modes, m = 0, ±2, ±4, . . . 
(lLongcope & Priesa l2007. and appendix). The disruption of 
the current sheet occurs at frequencies near and below the 
characteristic frequency, uja- For all values \m\ > 2 this 
evolution will be below the cut-off frequency so virtually 
no power will emerge from the inner region in any mode 
except the axisymmetric mode, m — 0. This means, as we 
demonstrate more precisely in an appendix, that the free 
energy in non-axisymmetric modes is directly dissipated 
by resistivity, and the axisymmetric energy alone leaves 
the X-point region. The axisymmetric mode will reflect off 
outer boundaries, i.e. the photosphere, to produce new waves 
propagating back into the X-point. These waves will not be 
axisymmetric, but will contain power in frequencies similar 
to those in the initial wave, eq. ( l24l . The non-axisymmetric 
components will therefore reflect from the X-point. We will 
henceforth focus on the axisymmetric, m — 0, in the inner 
region. 



The complex phase in the axisymmetric version expression 
of ( T32T ) represents the phase difference between the outward 
and inward waves, referred to position r = i n . Evaluation of 

that phase shows it to be nearly zero (i.e. J?(°) ~ e-™/ 2 "*) 
suggesting the waves reflect off the diffusive region. In order 
to provide coupling between inner and outer regions we intro- 
duce an axisymmetric reflection coefficient referred to radius 

n 

Ri(u;) = #°>(w)exp[2i(w/wA)ln(ri/4,)] • ( 34 > 

This will have the same magnitude, eq. d33l , but a phase rep- 
resenting the round-trip transit from fj to the diffusive region. 

An inward-propagating axisymmetric wave, with Fourier 
transform W-(uj), will produce an outward wave, W+ — 
RiW-. Using this in expression ( |23T l we find the total energy 
crossing the interface outward, 



P, dt 



1 

8u> A 



{Wl -Wl)dt 



1 

167T 



(i-|4| 2 )|^ 



dio 



(35) 



This negative value is the net energy dissipated at the X- 
point. Using eq. (l33l in this expression shows frequencies 
ui > uiaI^ to be damped most effectively at the X-point. 

3.4. Coupling to the outer solution 

The outer region, while not described by a simple X-point 
field, is nevertheless described by linear dynamics which can 
be summarized by a reflection coefficient, R , found in next 
section. An outward propagating wave, W+ (uj), will give rise 
to a reflected inward wave W- (lu); both waves are evaluated 
n. The outer reflection coefficient is then defined by 



at r 



W-(ui) 
W+{ui) 



(36) 



By direct analogy with Ri, we may associate 1 — \R \ 2 with 
the energy lost to the outer region, although here it will be 
by a failure to reflect (i.e. free radiation) rather than Ohmic 
dissipation. 

Coupling the inner and outer regions at the interface pro- 
duces a kind of cavity between the two reflectors. Requiring 
the waves in both expressions to match (i.e. the inward-wave 
from outside is equal to the outward wave from the X-point) 
gives the resonance condition 



Ro(u)Ri(u) 



1 



(37) 



This will be satisfied by a set of frequencies, u) n , constituting 
the resonant frequencies of the cavity. We have shown above 
that |i?i(cj)| < 1 for real, non-zero, frequencies; the exter- 
nal domain will satisfy the condition \R (ljM< 1 for real 
frequencies. The resonance condition, eq. d37l > will thus be 
satisfied only for complex frequencies. 

An example of this cavity analys is is provided by the analy- 
ses of ICraig & McClvmonj dl991l) and[Hassam| dl992l) . Their 
external potential field was a perfect X-point, just like the in- 
ner field. (In this simple case they did not need to separate the 
field into inner and outer regions.) The field was bounded by a 
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cylindrical boundary at r = L which was a perfect conductor: 
(L) = 0. This generates a perfect reflection 



R ((jj) — exp[ 2j(cj/u;a) ln(L/r,;) — i tt ] 



(38) 



where the — iir accounts for the inversion of U at the conduc- 
tor. Taking the phase of R^ ' to be zero, as it approximately 
is, gives a resonance condition for the cylindrical cavity 

Ro(u)Ri(u) =expj 2i(uj/uj A )[ln(L/e ri ) + in/4] -iir\ 
= 1 , (39) 
for Re (a;) > 00 The resonant frequencies are thus 



u>.„ 



41n 2 (L/^) 



(40) 



for n — 0, 1, The damping time of each mode is 

therefore proportional to the square of the logarithm of the 
Lundquist number, S = L 2 lua/v — L> 2 /£ 2 , as originally 
found by Hassam (1992). The real part of the resonant fre- 
quency is set by the requirement that an odd number of quar- 
ter waves fit between the node at r = L and the anti-node at 
r = £ v , separated by transit time \u(L/£ v )/uja- 
The quality of the cavity 

Re(ui n ) ihx(L/Q 2lnS 



Q 



Im(w„) 



(41) 



is roughly the number of times the wave must reflect off the X- 
point before it is damped. Lundquist numbers of astrophysical 
scale, S ~ 10 12 , would require a relatively large number of 
reflections off the perfectly conducting cylinder before the X- 
point would be able to dissipate all the energy Ohmically. 

4. THE DYNAMICS OF THE EXTERNAL FIELD 

4.1. Numerical solution 

We produce a numerical solution to the time-dependent 
equations of linear dynamics, eqs. ( fT3l and (fT4l . in the ex- 
ternal region. Diffusive effects are limited to the inner re- 
gion so we set 77 = in the external equations. The field 
Ai and U\ are represented on a uniform L y x L z rectilin- 
ear grid. Both fields are located on the same points, and the 
Alfven speed from the potential field, ua(x), is computed on 
these points as well. The Laplacian of A\ is computed us- 
ing a simple finite-difference operation which is second-order 
accurate on the uniform grid. Along the outer boundaries, 
y = ±L y /2, z = and z = L z , the conducting boundary 
condition U\ = is imposed. The grid is made large enough, 
compared to y\ and z x , that reflections from the upper and 
lateral boundaries do not return the X-point during a solution. 
The reflection from the z = boundary is the focus of the 
investigation. 

A propagating wave condition is applied at an interior 
boundary, r ~ fj. This is done by identifying two contigu- 
ous, closed paths of grid points called C_ and C+. Each 
path is approximately circular, at radius r_ and r + , such that 
(r+ + r_)/2 ~ Vi. The paths are nested, meaning no extra- 
neous grid points lie between them. Averages are computed 

3 This restriction follows from the desig nati on of inward and outward 
waves, as explained in the text following eq. (5TJ. 



along each path as a means of identifying the axisymmetric 
components and at both radii, 



tt5> . 



N± 



A+ = 



c = 



N± ^ 

± G± 



A l 



(42) 



where N± is the number of grid points along path C±. At 
each time step the values of A-, A + and U+ are computed 
from the present solution. These are used to set new values 
for U\ along the inner path, C_ . This is done using a finite- 
difference version of eq. ( f20b . 



-£/+-w A (V + +r_) 



A+ - A- 

r + — r_ 



+2W + (u) A t-1nr i ) 



(43) 

where W+ is the prescribed wave-form being fed into the so- 
lution. At each point along the inner path, C_, the value of 
U\ is set to the value U- computed for that step . The value of 
Ai along both paths is advanced using eq. ( fT4b , with 77 = 0, 
which requires only values of U\ at the same point — no in- 
formation from neighboring points. The solution on points 
inside the inner path is never used and is therefore irrelevant 
(the equations are solved there anyway since time-stepping is 
performed over the entire rectilinear grid). 

The procedure above imposes a prescribed, axisymmetric 
wave, W+, propagating away from the inner region and into 
the external solution. It also serves as a perfect absorber of 
axisymmetric waves incident on the inner region. The form 
of the incident wave passing out of the computational domain 
can be computed after the fact using the same finite-difference 
version of eq. d20i > 



W_(w A t-lttr<) - W++U-) - Wr++r_) 



A+ - A- 



r + — r_ 
(44) 

Since U\ is set to the same value at all points along C_, any 
wave component with to 7^ is perfectly reflected from the 
inner boundary. This is an approximation to the actual cou- 
pling, but as shown above it is a reasonably good one for low 
frequencies. 

The value of 77 appears nowhere in the equations of the outer 
region, and the solution is thus applicable to any value pro- 
vided £ v < rj. Its only effect on the outbound waveform, eq. 
d24l i. is in the delay time t . We eliminate this final depen- 
dance by setting 



W + (uAt -Inn) = 2u; A / C! 



1 



1 



ut A t 



t > 



(45) 

so the wave enters the external domain at the beginning of 
the solution, t = 0. When applying the numerical solution 
to a given 77, the current-sheet disruption will have occurred 
(i.e. 77 will have been switched on) at the negative time t = 
-ln(r|wA/?7)/2wA. 

The initial conditions are found by solving Laplace's equa- 
tion, V 2 Ai = numerically subject to the homogeneous 
outer boundary conditions and holding A\ at a constant value 
along the inner path C_ . This solution is then rescaled so that 



|(ri+rTj) 



Ax - A 
ri - r 



= -2 



(46) 



which gives a positive current of unit magnitude. The am- 
plitude of the linear solution is irrelevant, so in this way we 
normalize it to I cs . 
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Figure [3] shows a typical solution to the linear equations. 
This case uses a quadrupolar field with ^2/^1 = 4/5 and 
2/2/2/1 = —1/2. The source depth is set to d = O.lyi, and eqs. 
© and© give Zy — 0.4 yi and Bq — 1.92 t/>i/?/^. We use 
a 701 x 300 point grid with conducting boundaries at z = 0, 
z = 3yi = 7.5 z x and y — ±2.5 y\. (The grid spacing is 
therefore Ay = = 0.01 j/i.) The inner boundary is taken 
to be = 0.07yi; the inner path C_ is a white curve in the 
figure. The four panels show U%(y, z) at four successive times 
in the solution, depicting the initial departure of the wave (up- 
per left), its reflection from the photosphere (upper right), and 
then encountering the X-point upon reflection (lower panels). 
The central region is positive in the initial phases. 



CJ s t = 



a A t = 




FIG. 3. — The solution to linear dynamics amidst a quadrupolar potential 
field. Ui (y, z) at uipj, = 1.6, 2.2, 3.5 and 5.7, from upper left (color online). 
The white curve is the inner boundary C_ . Gray curves are the separatrices 
of the potential field. 

The inductive electric field Ui — vi ■ (x x B ) is the velocity 
perpendicular to qudrupolar potential field Bo- An axisym- 
m etric (m = 0) distribution of U\, such as the one predicted 
by Longcope & Priest ( 2007), corresponds to a quadruplar ve- 
locity field. Figure |4] shows the velocity vectors at two times 
in the solution. The plot assumes I cs > 0, so the initial cur- 
rent sheet was horizontal like the one shown in Fig. Q] (Had 
I cs < 0, the current sheet would be vertical, Ui < 0, and all 
the arrows would be reversed.) The quadrupolar velocity field 
is oriented in the same way as in steady-state reconnection: 
inward above and below the current sheet and outward from 
its ends. This sense eliminates the initial counter-clockwise 
azimuthal magnetic field perturbation Bi = WAi xx, thereby 
decreasing B = Bo + Bi in the inflow regions (they are rar- 
efactions) and increasing it in the outflow regions (they are 
compressions). Unlike steady reconnection models, the in- 
flow is not here imposed from the boundary; rather it arises 
from the rarefaction portion of the propagating FMW. 



4.2. WKB solution 

The structure of this typical solution can be understood in 
terms of an WKB approximation, valid for the components of 
highest frequency. For the spatially varying dispersion rela- 
tion, w(x, k) = ^a(x) |k , a ray path r(t) must satisfy Hami l- 
ton's equations (Kulsrud 2005; McLaughlin & Hood 2006a) 



dr 

di 
dk 
~dt 



duj 
dk 



k 

\k\ 



(47) 



(48) 



A ray is initialized with r(0) located somewhere on the circle 
r = Ti and k directed radially outward (its magnitude is ir- 
relevant). The full ray is then found by solving eqs. d47b and 
d48t . If r encounters the z = surface the ray is specularly 
reflected by changing (k y ,k z ) — > (k y ,—k z ) at that instant 
and the continuing the solution. The upper left panel of Fig. 
[5] show rays initiated at different positions along the interface 
circle. It also shows nested, closed wave fronts formed by all 
rays at the same time t. These approximate the leading edge 
of the U\(y, z) front, as for example in Fig.[3k. 
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FIG. 5. — Solutions to the WKB equations, eqs. J47M 48i. The separatrices 
of the potential magnetic field are plotted in each panel with light grey, (a) 
Rays and fronts for times t < tt r , before any reflection. The grey region 
includes the rays destined to reflect off the photosphere, (b) Fronts for times 
t > ttr. Broken lines show the sections reflected from the photosphere, solid 
shows the unreflected sections, (c) Rays (dashed) and fronts (solid) from the 
rays reflecting off the photosphere, (d) Rays (dashed) and fronts (solid) from 
the rays not reflecting off the photosphere. 

It is to be expected that the rays initially directed upward 
will not reflect from the photosphere, and this expectation 
is borne out. In fact only a fraction of those rays initially 
directed downward end up reflecting; these are within the 
shaded region of Fig. [5^. The remainder are refracted away 
from the large Alfven speed above the sources, and propagate 
upward in the end (Fig. |5j;). 

The first reflection occurs for the ray propagating straight 
downward, k = — z, after 



FIG. 4. — Velocity vectors corresponding to the electric field distribution, 
Ui(y, 2), at two early times in the solution. Arrows show the direction and 
magnitude of vi, although their length is capped at a moderate value for 
clarity. Dashed lines show the separatrices of the potential field and the solid 
curves are the wave fronts from the WKB approximation. 



*tr = 



dz 



1.42 



(49) 



(This is slightly less than the time, ln(z x /r^/wA — 1-74/u;a, 
which would be taken in a pure X-point field due the increased 
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field strength between the inner two sources.) The rays begin- 
ning with angles further from vertical reflect off the photo- 
sphere successively later (see Fig.[5jl), but only out to ~ 28° 
of downward (~ 15.5% of all rays). The ray bounding the 
reflected set is one that appears to graze the z — surface 
tangentially, and is therefore marginally reflected (k z — at 
the point of grazing, so its unimpeded motion is equivalent 
to specular reflection.) In this case the ray grazes a point be- 
tween the photospheric sources (compare its location to the 
separatrices shown in lighter lines). The ray is then refracted 
upward by the large Alfven speed of the outer source. 

The reflected rays are sent upward from z = 0, but only 
a fraction of these appear to focus back on the X-point. In 
the vicinity of the X-point, where vx(r) ~ WAf, Hamil- 
ton's equations are solved by rays forming logarithmic spirals 
(McLaughlin & Hood 2006a), 



it) 



Ti e 



LJ A(t — ti ) COS 9q 



fa[t) = <j)i+uj A (t-ti) sin 9 , 
(50) 

where 9o = tan (fc^/fc r ) is the angle of the ray relative to 
radial, and ti is the time it crosses r = at angle fa. Of 
the reflected rays, only those rays for which k r < 0, and thus 
cos #o < 0, will approach the X-point. The wave-front from 
these rays drapes around the X-point, later wrapping about it 
as the rays complete spiral circuits. Initial draping is evident 
in Fig. |5p. 

The set of non-reflected rays, bounded by the grazing rays, 
forms a locus of rays bound upward. This WKB construction 
is reflected in the structure of the full solution at late times, 
as shown in Fig. [6] A concentrated pulse of inductive elec- 
tric field, U\, approximately coincides with the non-reflected 
rays. The wave-front of the reflected rays adjacent to the graz- 
ing ray (dashed curve), forms an approximate lower bound to 
this region. At this stage the inward flow has become decou- 
pled from the X-point, so it is no less clearly related to the 
reconnection inflow. The current density, J x = — V 2 Ai/47r, 
shown in Fig.|6j3, is also concentrated within this region. Be- 
low the reflected front the current density is negative, while 
above it the current has the same sign as the current carried 
by the sheet: J cs > 0. The total current within the positive 
concentration slightly exceeds, that from the initial sheet. 




FIG. 6. — The solution at a time later than those shown in Fig. [5] The left 
panel shows U\ (y, z) in an inverse grey scale, but expanded by a factor of 
3. Superimposed on this are the wavefront, from the same time,of all non- 
reflected rays (horizontal solid segment) as well as the grazing rays (vertical 
solid segments). Dashed segments are a small section of the wave front from 
reflected waves. The right panel shows the current density J z from the same 
time, with the same superimposed WKB solution. 



4.3. Reflected waves 

The WKB solution shows that only a small fraction of ray 
paths (15.5% in this case) leaving the X-point vicinity reflect 
from the photosphere. Of these an even smaller fraction con- 
verge back on the X-point. All but one of these rays approach 



the X-point in a spiral, and therefore approach more slowly 
than they left. We therefore expect the reflected wave inci- 
dent on the X-point, W- , to be of lower frequency than that 
emitted, W+. The WKB approximation is, however, applica- 
ble only to very high frequencies. We therefore look to the 
full numerical solution to understand the waves incident on 
the inner radius . 

Figure|7]shows the incident wave, computed from eq. d44l ). 
in the same solution from Figs. [3] and [6] It begins to depart 
from the initial value, W- = 2uj a I cs , at the time 2t tI , taken 
for a vertical ray to make a complete round-trip. Its evolution 
following that is on time scales longer than the initiating wave 
( W+ shown dashed), as predicted by the WKB solution. 




FIG. 7. — The waves at the interface from the solution shown in Figs. [5] 
and [6] The lowest panel shows W— (solid) and (dashed), normalized 
to 2u)^I cs . Diamonds show the times of the the four panels in Fig. [5] and 
a square shows the time of Fig. \§\ Vertical dotted lines show the time of 
initiation, t = 0, and first reflection, t = 2tt T - The middle panel shows 
the outward energy flux crossing the r = rj boundary, computed using eq. 
423t and scaled to The dashed curve shows the version for purely 

outging wave, plotted in Fig. [2] The top panel shows the inner current, /; 
(broken), and flux function, A\ (solid), both scaled to their values at r = r^. 
A dashed curve shows the current for the purely-outgoing case. 

While the high-frequency components are only partially re- 
flected, it seems the lowest frequency components are more 
completely reflected: W- — > —W + by the end of the solu- 
tion. This was predictable from the outset. In order to sustain 
the initial equilibrium, where W— = — W+, it is necessary for 
the photosphere to be a perfect inverting reflector at cj = 0, so 

i? o (0) = — 1. Nor can this result from the outer computational 
boundaries, since a current sheet equilibrium is achievable, in 
principle, in an unbounded half-space. 

As a consequence of the high reflectivity at low frequen- 
cies, the Poynting flux through the interface (middle panel) 
approaches zero. The net current inside the inner region, 



T - 1 dAl 

1 ~ 2^ dr 



1 



4wa 



(W_ - W+) 



(51) 



reverse sign and achieves a magnitude close to its initial: /j — > 
— J cs . This same behavior was observed in the case of a per- 
fect cylindrical reflector, where it was fou nd to eliminate and 
then r everse the residual X-point current ( Longcope & Priest! 
120071) . It would establish a new vertical current sheet in place 
of the old horizontal one, if there were not resistivity in the 
X-point vicinity. The reversed current also leads to a reversal 
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of the flux function, Ai(rj). This tracks the current as pre- 
dicted by equilibrium eq. dSj, except for the dynamical effects 
of waves crossing the interface. 

Based on the current Ii it would seem that the external sys- 
tem has been returned to its initial state; this is the case for 
the perfect reflection from a conducting cylinder. The free 
magnetic energy and kinetic energy of the external solution 
are 



E M = ^~ [ \VA 1 \ 2 dydz , 
8tt J 



1 f U 2 
Efc = — I — dy dz 



8tt 



(52) 
(53) 



T~>ri 



The energies from the foregoing solution are plotted in Fig. 
[8] along with the total energy, Em + Ek (dashed). The ini- 
tial magnetic energy equals the value from eq. ([Tol l, and all 
energies are scaled to this. Since resistivity is dropped from 
the external dynamics, the total energy changes only by en- 
ergy flux, Pi, across the inner boundary. (Since U\ = on all 
outer boundaries there is no Poynting flux across them.) 
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FIG. 8. — The energetics of the external dynamics. Solid curves in the bot- 
tom panel shows the kinetic (lower) and magnetic (top) energies, scaled to 

AB^° ut '. Diamonds show the times of the the four panels in Fig. [5] and a 
square shows the time of Fig. [6] The dashed curve is the total energy. The 
middle panel shows the Poynting flux (solid, scaled to 2uj\I^ b ) and its inte- 
gral (dashed, scaled to Ai?j^ ut '). The top panel shows the same quantities 
computed for the case of perfect reflection from a concentric cylindrical con- 
ductor. 

As the FMW passes into the exterior region, it draws en- 
ergy inward, driving the total external energy down. The 
kinetic energy increases, but the magnetic energy decreases 
more rapidly. The incidence of the reflected wave, W-, on 
the interface reverses the sense of energy flux, since it causes 
\W—\ to decrease. At the instant |W_| — |W+|, the sense of 
energy flow across the interface changes from inward to out- 
ward. Owing to the different time scales of the outward wave 
and reflected wave, the positive energy flux lasts longer than 
the negative. The result is that the total energy of the external 
solution increases by ~ 20% of its initial value, by the time 
the reflected wave is fully absorbed. 

The energy increase contrasts with the cases of perfect re- 
flections, previously studied (shown in the top panel), where 
the outgoing and reflected time-scales match and the net en- 



ergy flux approaches zero. In the more realistic case consid- 
ered here the state following absorption of the reflected wave 
includes an upward-propagating FMW (see Fig. [6]) whose ki- 
netic energy remains in the solution. At the time of Fig. [6j 
the kinetic energy density, U 2 /v\, is evidently concentrated 
within the upward pulse confined to the WKB envelope. It 
seems that the slight reduction in equilibrium magnetic en- 
ergy accompanying the slight reduction in current < I cs ), 
is roughly compensated by the magnetic perturbation of the 
FMW. 

This corroborates our initial expectation that imperfect re- 
flection from the photosphere would reduce the ability of the 
X-point to directly dissipate energy. After a single reflection, 

the FMW has extracted ~ 0.2 AE M ut ^ from the interior cur- 
rent sheet. This is energy not available for direct dissipation, 
later, at the current sheet. It is the first part of the total energy 
which will be radiated by fast magnetic reconnection. To es- 
timate the effect of repeated reflections, and thereby estimate 
the total energy radiated in FMWs, we turn to the reflection 
coefficient R , defined in eq. (|36"1 ). 

5. LONG-TIME SOLUTION: REPEATED REFLECTIONS 

5.1. Repeated reflections 

The analysis above showed how the energetics of reconnec- 
tion could be computed using only the amplitudes of incom- 
ing and outgoing waves at the interface r — fj. These waves 
satisfy linear equations so they can be computed from a su- 
perposition of successive reflections 



W ± = wf 



w. 



(1) 



W. 



(2) 



(54) 



The leading terms, W± , are the solutions computed in the 

previous section: wi. is the wave from the current disruption 

given by d45l l. and and is its reflection from the external 

solution shown in Fig. [7] The initial wave form, , is valid 
only when the incident wave is steady: W- = 2I CS . We thus 

define the next correction, W_l , to be the response to the 

departure of the incident wave from this, W { _ 0) - 2/ cs . The 

Fourier transform of this departure is — AttI cs S(uj), and 
its reflection from the X-point is 



(w) = Ri (w) [ W w (w) - 4irl cs 5(u) 



(0), 



(55) 



after using the fact that Ri(0) = 1. This correction wave 
then reflects from the exterior to provide the correction to the 

incident wave W_ — R W^\ The remainder of the terms 
are subsequent reflections of the correction 



(n) 



= R Wi n) 



n > 2 . (56) 



The X-point reflection coefficient is given in analytic form 
by eqs. d32b and 041 ). Its dependance on 77, through the phase- 
delay factor in eq. ( f34b . is the only place where resistivity 
enters the computation. The coefficient of outer reflection, 
R (uj), can be found from the numerical solution of a single 
reflection, such as the one from the previous section, by tak- 
ing the ratio of Fourier transforms W+ and W- . It proves 
to be more efficient to perform the computation with a dif- 
ferent input wave form, W+, whose power is more evenly 



11 



distributed over frequencies. For this we use a Gaussian 
W+(wa£ — lnr*i) = exp[— t 2 /t 2 ]. Taking r somewhat larger 
than 1 /u)a samples a broader range of frequencies, and more 
equitably, than the smoothed-out step of eq. (l45l l. We solve 
for the linear response as before, except the initial condition 
is A x = 0. 

Figure [9] shows the magnitude of the reflection coefficient, 
\R a \, computed this way. The reflected wave had decayed 
to nearly zero, and the run was stopped, before any effect 
had reflected from the upper (z — L z ) or lateral boundaries 
(y = ±L y /2), so the response is applicable to arbitrarily large 

domains. |i? D | is slightly larger than the coefficient of the X- 
point, \Ri\ (dashed), given in eq. d33l l. It is largest at fre- 
quencies below that from the transit time, ui = 7r/2i tr , and 
approaches unity at cj —> 0, as it must. At higher frequencies, 
R approaches a floor near the level of reflection predicted by 
the WKB calculation: 15.5%. 



structure of | R\ | and \R \. Were it not for the broadening the 
composite wave (i.e. eq. l54l ) would be a series of steps. In- 
stead it becomes increasingly sinusoidal as higher frequency 
components are damp ed away at each reflection. (The same 
tendency was found by Longcope & Priest! ((2007), and is evi- 
dent in their Fig. 11.) 
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FIG. 9. — The external response to a Gaussian wave. The top panel shows 
the input, W+, a Gaussian (dashed) and the response W— (solid). The ver- 
tical dotted lines show t = and t = 2tt r - The bottom panel shows the 
magnitude of the reflection coefficient |JJ | (solid), and the analytic form of 

the X-point reflection, \Ri | (dashed). The vertical dotted line is ui = 7r/2tt r , 
and the horizontal dotted line is at 15.5%, the fraction of rays reflected in the 
WKB solution. 

The X-point reflection, Ri, has a finite negative slope to 
the right of uj = 0, and by realizability conditions, a positive 
slope to its left. Its first derivative is discontinuous at ui = 0. 
Such a discontinuity is accompanied by an impulse response 
decaying as t" 1 . This residual response was first noted by 
Hassam (1992), and is related to asymptotics of eq. d45l l due 
to the shrinking region of residual current ( Longcope & Priest! 
120071) . The external reflection coefficient, by contrast, appears 
to approach w = with zero slope. The assymptotics of the 
external region will be governed by reflection, with losses, 
from the photosphere. 

Figure [10] illustrates this procedure for a Lunquist num- 
ber S = z 2 x /£ 2 = 3 x 10 10 . This value enters only 
into the phase factor in eq. d34i i. representing the transit de- 
lay, ln(j-j /£,,)/£*; a = 10/wa, between interface and resistive 
scale. This plus t tI is half the round-trip delay for reflec- 
tions off both the photosphere and the X-point. It is there- 
fore approximately half the delay between the rising edges 

of W± ^ and W±\ Each reflection is a step of opposing 
sign, delayed by this amount and broadened by the low-pass 



20 40 60 SO 100 120 140 20 40 60 SO 100 120 140 

FIG. 10. — The different components of the reflecting waves, for Lundquist 
number S = 3 X 10 10 . Left and right panels show W+ and W— respectively. 

(n) 

Dashed curves are individual components, W± for n = 0, 1, 2, and 3 (from 
bottom to top). The solid curve is the sum up to n = 14. The top axis shows 
time in units of the full transit, tti + ln(rj/£n)/wA. 

The energy flux due to this wave, plotted in Fig. [TTJ shows 
the initial behavior found in Fig. [8] a brief positive flux, fol- 
lowed by a longer negative (inward) flux. This repeats once 
more when the wave reflects from the X-point, this time with 
a longer negative portion and shorter positive portion, lead- 
ing to less net outward flux. Thereafter the waves become 
increasingly sinusoidal assuming a phase relation where the 
upward zero crossing of W+ is slightly ahead of the down- 
ward crossing of W- \ the ratio is ~ £tr w A/ ^(ri/£ v ), due to 
the difference in transits. This phasing leads to a steady de- 
crease in the net energy transport until it converges on a value 

~ — 0.41A£^° ut '. Free magnetic energy from the external 
field is, therefore, moved into the inner region and dissipated. 

The remaining free energy, ~ 0.59A£<° ut) , is therefore con- 
verted into FMWs propagating away from the reconnection. 
The full solution (i.e. Fig. [3]), covering only the first oscilla- 
tion, showed a single wave accounting for about a third of the 
total. It seems that new waves are launched under repeated re- 
flections, each with smaller amplitude. The total energy in all 

waves is the - 0.59A£4° ut) carried away from the exterior. 
5.2. Scaling with Lundquist number and free energy 

The external solution, leading to R (uS), does not depend 
on the resistivity. The inner reflection, Ri, involves resistivity, 
through^ = y/r)/uj A , in the phase-factor e H"/"A)Hn/tr,) 
in eq. (l34l >. It is therefore easy to perform the same cal- 
culation for any value of r\. The results all resemble Fig. 
[TTI with predicable differences. As rj decreases the transit 
time, ln(ri / £ v ), increases and the oscillation period increases. 

Successive reflections, W^™ -1 " 1 and W±, are spaced farther 
apart therefore cancel each other less resulting in more grad- 
ually decreasing total waves, W±. The phase difference be- 
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FIG. 1 1. — The energy flux across the interface from the waves in Fig. 
|10| The bottom panel shows W-y. (dashed) and W- (solid); this extends, 
through independent computation, the bottom panel of Fig. [7] The middle 
panel shows the Poynting flux, Pj and its integral, AE in the same format as 
in Fig.[8] The top plot is current Ii. 

tween these approaches 180° since ttr becomes a diminishing 
fraction of the whole period with smaller r\. This results in a 
diminishing net transport inward and thus a smaller peak val- 
ues of Pi. This lower mean power persists for a longer time 
before the component waves decay. The result is a net energy 
transfer which depends only very weakly on the logarithm of 

S: from -0.50A£<° ut) at S = 3 x 10 5 to -0.35A^° ut) at 
S = 3 x 10 12 . A sweep of this range is well fit by a fraction 
-0.6 + 0.0076 In S. 

We therefore find, neglecting the extremely weak depen- 
dence on S, that a fixed amount of magnetic energy is con- 
verted to FMWs by the reconnection. That energy is less than 
all the free magnetic energy outside the radius r, = 0.0 7yi we 
used for our computation. Using eq. ( [Tol l, we find the wave 
energy to be approximately equal to all of the free magnetic 
energy outside the larger radius ~ 0.19j/i. This radius is 
outside the validity of our calculation since Ao(y, z) inside 
such a large r, cannot be approximated by an X-point, as in 
eq We extrapolate to these values in order to illustrate the 
expected behavior of the system with significant free energy. 
In such a case, illustrated in Fig. [12] all of the free magnetic 
energy inside the region would be dissipated by the resistivity, 
all the energy outside would bet carried away by FMWs. 

This foregoing scenario appears to be virtually the same 
regardless of the value of resistivity responsible for the dis- 
sipation, provided it is sufficient to cause fast disruption of 
the current sheet. The total free energy in the initial field de- 
pends on the length of the current sheet, 2A (i.e. through eq. 
13), which itself depends on the current I cs . The net cur- 
rent in the sheet is in turn set by the flux discrepancy, Aip, 
imposed by ideal evolution driven from the boundaries (for 
instance through flux emergence). The right side of Fig. Q~2] 
shows the total free energy (solid line) and the amount con- 
verted to FMW; the upper panel gives the fraction of total 
energy converted to waves. That fraction becomes smaller as 
the initial energy becomes small. Those cases represent sheets 
with very small currents which are therefore themselves phys- 
ically small (A <C y\). The free energy becomes increasingly 
concentrated in the small neighborhood of the current sheet, 
where it is more easily dissipated. So dissipation eliminates a 
larger fraction of the energy, which is itself smaller due to the 




FIG. 12. — Left: The quadrupolar field with y^ /y\ = —1/2 and ^2/^1 = 
4/5, used in the forgoing calculations. The dashed circle in the bottom panel, 
and the grey area in the top, define the interior region whose energy is entirely 
dissipated. The energy initially outside is converted to FMWs. The bottom 
panel shows the field lines, contours of A(y, z), and the top shows lines 
of the non-potential field, contours of A±(y, z). Right: the free magnetic 
energy (dashed) and energy carried by waves (solid) versus the change in 
flux Aip/ipi. Plotted is the logarithm (base 10) of the energy in units of 
ipf . The top panel shows, as a fraction of the total free energy, the energy 
ultimately converted to waves (solid) and the energy initially present in modes 
with m > 2 (dashed) and then directly dissipated. The axis along the top is 
the logarithm of the maximum permissible Lunduist number, from minimum 
permissible resistivity, r)t = lo\A . Diamonds correspond to the current 
sheet shown in the right panel. 

smaller current. 

The foregoing analysis considered only axisymmetric 
waves, m — 0, in the vicinity of the X-point. The initial cur- 
rent sheet also included modes with m > 2, but it is shown 
in the appendix that the energy in these modes is dissipated 
immediately by the resistive diffusion. The energy in these 
modes, ~ 0.27 I 2 S , is included in AEm, of eq. ©, and pro- 
vides a lower limit on the energy directly dissipated by Ohmic 
diffusion in our model. Plotted as a dashed curve in the upper 
right panel of Fig. [12] it accounts for between 3% and 11% of 
the total. 

Releasing significant free energy requires that the initial 
current sheet have non-trivial physical extent. In this case 
our uniform resistivity, 77, must be interpreted as the result 
of complex (turbulent) small-scale physics occurring rapidly 
within the current sheet. In order to simplify our analysis 
we have assumed this turbulent diffusion is uniform in space 
and constant in time, following the onset of the turbulence. It 
must also be large enough to disrupt the entire current sheet 
over a time comparable to the Alfven transit across the sheet: 
r\t ^> waA 2 , which is equivalent to our earlier assumption 
that £ v ^> A. This requirement places an upper limit on the 
turbulent Lunquist number — a limit which depends on I cs 
through A. The logarithm of this limiting value is plotted 
across the top right of Fig. [T2] It is evidently necessary to 
have S ~ 100 in order to rapidly reconnect one percent of the 
photospheric flux; a current sheet like that in the left panel of 
Fig-El In this case 40% of the free energy released would be 
converted into FMW, and the remainder would be dissipated 
by the turbulent processes at the X-point. 

6. DISCUSSION 

We have used a simplified model to study the response of 
the large-scale coronal magnetic field to reconnection at a cur- 
rent sheet. The reconnection reduces (almost to zero) the cur- 
rent in the sheet by transferring magn etic flux across it. As 
reported by Longcope & PriesJ (120071) . the current removed 
from the sheet is carried away at the front of a fast magne- 
tosonic wave. A portion of this wave reflects from the pho- 
tospheric boundary, and is refracted back toward the X-point. 
Subsequent reflections of the wave between the X-point and 
the photosphere lead eventually to the elimination of all cur- 
rent, and with it all free magnetic energy. 
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The number of reflections required to reach the potential 
field is roughly the logarithm o f the global Lundquis t num- 
ber, corroborating t he findings of Craig & M cClvmoml d 199 lb 
and Hassan] (119921) . The total reconnection time is therefore 
~ In 2 S times a single Alfven transit (i.e. 1/cja)- This final 
approach resembles the case of perfectly reflecting concentric 
cylindrical boundary because the photospheric boundary has 
a much higher reflection coefficient than the resistive X-point 
at very low frequencies. During early reflections, however, 
the FMW contains frequencies above uja/^S, of which a 
sizable fraction are directed vertically upward from the photo- 
sphere. We thus find, at least for this simplified model, that an 
appreciable fraction of the initial free energy is carried from 
the current by fast magnetosonic waves; in the case we ex- 
plored in detail 25% - 60% of the energy was converted to 
waves, depending on the size of the initial current sheet. 

We expect some of the basic elements of this result to hold 
in models more sophisticated that the one we used for our 
detailed study. A common element in all models of fast mag- 
netic reconnection is that the reconnection electric field, i.e. 
the flux transfer, occurs on small scales, while free mag- 
netic energy is stored over vastly larger scales. Information 
about the flux transfer, and associated current reduction, must 
be transmitted to the larger corona. This includes signifi- 
cant volumes which are not magnetically linked to the re- 
connection site, for which transmission must occur through 
fast magnetosonic modes. Previous studies of unsteady fast 
reconnection have shown fast magnetosonic waves propagat- 
ing into the unreconnected flux to be responsible for creating 
the reconnection inflow, the same inflow ass umed as a bound- 
ary condition or driver in steady models dLin & Led Il994t 
iHevn & Semenovlll996t iNitta et alj|200lb . Our model con- 
tains similar fast mode rarefaction waves, and reveals the frac- 
tion of energy they contain. 

One of the most dubious simplifications is our use of uni- 
form classical resistivity, 77, as the means of generating a re- 
connection electric field. The simple mathematical form of 
this effect permits a more thorough analysis than would more 
complex physical effects. Since our goal was to study the re- 
sponse of the large-scale field, where the electric field is irrel- 
evant, we chose the simplest possible form for the small-scale 
effect. We expect that a more accurate treatment of the small 
scales would reveal complex flows, whose overall effect might 
be characterized as an effective diffusivity. It is this kind of 
turbulent 77 we use in applying our results to the reconnection 
of global currents, as in Fig. [12] Were the turbulent diffusivity 
computed with any degree of self-consistency, it is unlikely to 
remain constant through the many repeated reflections of the 
FMW. Remarkably, we find that its actual value has no effect 
on the magnitude of the X-point reflection (see eq. l34l ). crit- 
ical to the ultimate energy dissipation. Thus we expect the 
basic scenario revealed in our simple model would be found 
in more sophisticated ones. 

The equations governing the large-scale response, namely 
eqs. ( fTTT i and ( fT2l . are linear, purely two-dimensional (with no 
magnetic field in the ignorable direction; no "guide field") and 
assume zero pressure. The assumption of linearity is the most 
easily justified, since the initial current sheet contributes a 
small c orrection to the field far from itself. Longcope & Priest! 
(2007) analyze in detail the requirements for linearity in the 
internal solution, and find that it is a reasonable approxima- 
tion for many times 1 /wa, by which time the reflected wave 
comes into play. 



Significant pressure can alter the nature of the equilib- 
rium from which the energy must be released. A recent 
numerical solution found that only a tiny fraction (~ 3 x 
10 -9 ) of the magnetic energy in a pressure-dominated equi- 
librium was released by sudd enly enhancing the resistivity 
(iFuentes-Ferandez et alJ l2012). They found that resistive dif- 
fusion does redistribute the current, as in the /3 = case, 
but that the new distribution of Lorentz forces are compen- 
sated by a new distribution of pressure. This redistribution 
process is localized to the current sheet so that instead of an 
axisymmetric FMW reducing the free energy, a small wave 
with m = 4 is launched, with negligible energetic con- 
sequence. A still more recent experiment at lower values 
of f3 continued to observe reconnection leading to pressure 
re-distribution rather than a de crease in the sheet's current 
( Fue ntes-Fernandez et al.ll2012l) 

At th e other end of the process, iMcLaughlin & Hoodl 
(2006b), found that pressure in the vicinity of the null kept 
the fast magnetosonic speed above zero, thus mitigating the 
focusing of wave energy there. Pressure thus reduces the 
amount of energy initially released and the fraction of that dis- 
sipated. By focusing on the zero-pressure limit we are finding 
the maximum possible energetic effect in solar flares, which 
pressure would almost certainly reduce. 

We believe our assumption of two-dimensionality places 
the most severe limitations on the applicability of our model. 
In the absence of a guide field (B x ) all shear Alfven modes are 
polarized in the ignorable direction (x) and are completely un- 
coupled from the reconnection dynamics. It is for this reason 
that the response propagates isotropically from the X-point as 
a FMW. Adding any amount of an ignorable component to the 
system will result in a global response including both FMWs, 
propagating in all directions, and shear Alfven waves con- 
fined to the separatrix field lines. While the FMW will reflect 
from the photosphere specularly, as they do in our model, the 
reflected Alfven waves will travel back along the same field 
lines to reach the X-point again. This far more complex, and 
more interesting, system will need to be studied in the future. 
The present study can be taken as a limiting case, and can pro- 
vide an outline, for such a future investigation. We expect the 
FMW component of that system to behave in a manner sim- 
ilar to our model, but to account for a smaller portion of the 
energy, since Alfven waves will contain some as well. 

Many past observations have provided evidence for both 
FMW and shear Alfven modes initiated by solar flares. Ev- 
idence of fast mode shocks are fo und in coronal type II ra- 
dio bursts (Pavne-Scott et al. 1947) an d chromospheric More- 
ton waves (Moreton & Ramseyl[l9 60: Moreton 1960), while 
shear Alfven waves are observed in post-flare loop o scilla- 
tions dAschwanden et al.|[l999t [Nakariakov et all 119991) . In 
both cases their timing strongly suggests these disturbances 
are initiated by flares or by CMEs, although the details of the 
initiation remain unclear. The wave front in our model in- 
cludes a "skirt" along the lower boundary , evident in F ig. [5b, 
which might manifest as a Moreton wave (Uchida 1968J). The 
upper portion of the wave, if correctly oriented, could shock 
to produce metric type II radio signatures. 

A common model of these observed waves has been that 
direct energy dissipation creates a p ressure pulse driving a 
FMW outward in all directions (see Vrsn ak & Cliverl 120081 
and references therein). Our model provides a detailed picture 
of the initiation mechanism which differs in several important 
respects form pressure-pulse models. Our primary conclusion 
is that, contrary to the common assumption, the broad initial 
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distribution of free magnetic energy makes it impossible to 
create a local pressure pulse through its rapid dissipation. We 
find instead that rapid diffusion results in a current redistri- 
bution whose newly unbalanced Lorentz forces produce the 
FMW. We admittedly discard the pressure which would make 
possible the alternative wave-generation mechanism. Were it 
included, however, the thermal energy could not exceed the 
magnetic energy directly dissipated. We find this to be 3%- 
1 1% of the total, and therefore smaller than the wave energy 
produced through Lorentz forces. 

The plasma flow direction is one potentially observable dif- 
ference between our model and the pressure-pulse waves pro- 
posed previously. The pressure pulse or piston driver leads to 
outward (compressive) flows in all directions, while our wave 
has both outward and inward (rarefaction) flows along differ- 
ent portions of the cylindrical front. The outward portion is 
an extension of the inner reconnection outflow jet and thus 
related to the fast m agnetosonic termination shock predicted 
in various models (Forbes 1986). The inward portion is the 
reconnection inflow, but its speed is a function of the external 
field rather that the inner reconnection rate. The relative lo- 
cations of these two components depend on the orientation of 
the initial current sheet, and thus on the sign of / cs . 

Moreover, since all its energy is introduced at the outset, 
the pressure-driven wave has an amplitude that decreases with 
distance from the source. In contrast, the FMW continues to 
draw energy from the large scale magnetic field and thus de- 
cays much more slowly; in the vicinity of the X-point its am- 
plitude rem ains constant even as its n et energy increases ex- 
ponentially (Longcope & Priest 2007). It remains to compare 
the amplitude profiles of observed FMWs to the predictions 
of these different models. 

Finally, our model shows how the dissipation of magnetic 
energy at the smallest scales, i.e. the current sheet itself, re- 



quires repeated reflection from distant boundaries (i.e. the 
photosphere), and is therefore not entirely local. The initial 
phase of flux transfer dissipate s only a very s mall f raction of 
the energy, as found by Longcop e & Priest! (120071) . and ex- 
pected from arguments based on the finite energy density and 
very small volume. Reflected waves will focus back onto the 
X- point due to its vanishing Al fven speed. As firs t predicted 
bv lCraig & McClvmontl d 199 lb and lHassaml dl992l) . this fo- 
cusing leads in the end to the dissipation of significant energy 
within the small region after many reflections. The energy 
dissipation thus persists far longer than the time taken for the 
initial flux transfer that we call reconnection. 

It is tempting to see in this persistent dissipation a possible 
explanation for flare durations longer th an free coolin g fol- 
lowing a single impulsive energy release dWarrenl 2006). Be- 
fore doing so we must confront the effect, alluded to above, 
of three-dimensional geometry. With a guide field component 
the fast magnetosonic speed no longer vanishes at the X-point, 
and we would expect far less of the wave energy to focus back 
there. Alfven waves, on the other hand, will be confined to 
the field lines and will thus reflect repeatedly back to the re- 
connection site. This kind of wave reflection and dissipation, 
sometimes invoked to explain the persistent energy release in 
flares, is indirectly related to the fast mode version we have 
studied. Genuine fast mode focusing would, however, occur 
at a coronal null point in a three-dimensional magnetic field, 
as it does in our two-dimensional version. 
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APPENDIX 

NON-AXIS YMMETRIC WAVE COMPONENTS FROM CURRENT SHEET DISRUPTION 

The flux function for the straight current sheet of length 2 A, shown in Fig. [T3h . can be expressed in pol ar coordinates as an 
expansion valid either for radii r < A or for radii r > A. The external expansion, for r > A, was given in (Longcope & Priest! 
I2007FI 

evens , \\\\ 

A{r,<t>) = \B>^ cos(20) - 2/ cs mM 2 ^) + 4/ cs V 7 =- V) " f^f cos(m0) , (Al) 
2 V A/ ' m(m + 2)!! V r J 

where (m + 2)!! = 2 • 4 • • • m ■ (to + 2) and (to — 1)!! = 1 • 3 • • • (m — 1). The first term in ( IA1 b reproduces the simple X-point, 
while the remaining terms form the perturbation, Ai(r, <f>), due to the current sheet, shown in Fig. [13b. Expression ([8j includes 
only the leading term in the perturbation term (axisymmetric: m = 0), while we consider here the contributions of the remaining 
terms, m > 2. It is evident from Fig. [13b that this is dominated by m = 2, but also includes appreciable contributions from 
m = 4, 6, and so on. 
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FIG. 13. — The field in the vicinity of the current sheet, rendered by contous of the flux function A(y, z). (a) The flux fun ction A(y, z), (b) the pertu rbati on 
Ai(y, z) and (c) the non-axisymmetric contribution from m > 2. The dashed circle at r = A divides the external expansion )Alt from the internal one {A3) 

The internal expansion (r < A) of the entire flux function is 

evens 

A(r, <j>) = ±,4 (0) (r) + A(m) M C0S (W>) ( A2 ) 

m=2 

where m is once again even and 

odds 

A(m) M _ 16/cs y nil (r\ n 

7T j^j (n 2 -2n)(n 2 -m 2 )(n-l)!!VA/ ' ^ ; 

(by convention we take 0!! = 1). While the external series involves only a single inverse even power of r for each mode, the 
internal expansion involves all odd powers for each mode number m. These decrease asymptotically as n~ 4 and the series 
therefore converges rapidly. From these expressions we find the total free energy contributed by non-axisymmetric modes (to > 
2), the field in Fig. [T3b to be, 0.27/ 2 s . This contribution is included in the total energy, eq. (0, but is a very small fraction 
provided A -C z x . 

Starting with the linearized resistive MHD equations (flTT) and (fT2l . we again take a polar expansion of our perturbed flux 
function and electric field variable to obtain the scalar eqs. ( fT3l ) and (14\ . only now keeping the resistive term: 

9U±^^lr^ Al (A4) 

9A = -U 1+V V 2 A 1 . (A5) 
at 

Introducing the operator B m = (r + m) allows us to write the laplacian as V 2 = r~ 2 B_ m B m . We may then define the 
scalar function C^ n) = B m yl( m ), and express the azimuthal terms of eq. (IA4l i as 

^-U^ = -ujiB_ m C^ (A6) 

Operating with D m on eq. ( IA5I ). gives 

= -B m U^ + r?B m (^D_ m cW) (A7) 
Changing variables to s = ln(r/A), and letting /' denote df /ds, we arrive at these telegrapher-like equations: 

— C(" 1 ) = - m C/( m ) + ^e" 2s C" (m) - 2C" (m) - (to 2 - 2to)C ( " i) (A8) 

!L U {m) =-lo\ (C" (m) - mC (m) ) . (A9) 

The uniform resistivity 77 is tu rned on at t — 0, disrupting the current s heet . We numerically solve the equations of motion using a 
finite difference method: eq. (IA9I > is updated explicitly, while for eq. (I A8t we employ an operator splitting method, updating the 
first two terms explicitly and the diffusive term implicitly. We have a stationary initial condition, where U <<rn \s, 0) = 0. C^ m \ 
derived from the perturbation Ax, has r < A initial condition 

00 

aM(Sj Q) = 16ics g ^ !! e « s , (A10) 

7r (n z — 2n){n — m)(n — 1)!! 



while outside the current sheet 



C^ m \s,0) = \ / cs ' Z>2, (A11) 



since B m r~ m = 0, for m > 2. This regime well reproduces Fig 3. from Longcope & Priest (2007) for the m — mode. The 
to > 2 modes, however, are rapidly driven towards zero, even when 77 = waA 2 , as seen for m = 2, 4 in Fig. [14] The modes 
do propagate outwards for longer times, though not at a discernible level. Larger to modes are driven towards zero progressively 
faster. Magnetic energy contributed by the terms to > 2 is therefore dissipated locally by the resistivity, while the axisymmetric 
contribution leads to FMWs. 
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